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Abstract: We consider the Newtonian dynamics of a massive particle in a one dimem- 
sional random potential which is a Brownian motion in space. This is the zero temper- 
ature nondamped Sinai model. As there is no dissipation the particle oscillates between 
two turning points where its kinetic energy becomes zero. The period of oscillation is a 
random variable fluctuating from sample to sample of the random potential. We com- 
pute the probability distribution of this period exactly and show that it has a power 
law tail for large period, P{T) ~ T~^/^ and an essential singluarity P{T) ~ exp(— 1/T) 
as T ^ 0. Our exact results are confirmed by numerical simulations and also via a 
simple scaling argument. 
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Transport in random systems is an important and common physical phenomenon. 
Even small amounts of disorder in physical systems can dramatically modify their static 
and dynamical properties and lead to interesting and complex new physics. In presence 
of quenched, or time independent, disorder it is necessary to average physical quantities 
over different realizations of the disorder. This disorder averaging is often technically 
very difficult to carry out and few exact results exist, even in one dimension. In this 
Letter we study the zero temperature frictionless Newtonian dynamics of a particle in 
a one dimensional random potential which is a Brownian motion in space. Due to the 
absence of friction, the particle oscillates back and forth between two turning points 
where the kinetic energy becomes zero. The period of this oscillation is a random 
variable, varying from sample to sample of the disorder. The main result of this Letter 
is to present an exact result for the distribution of this time period. 

In general the dynamics of a particle in a one dimsional force field can be represented 
by the Langevin equation 



where 7 is the friction coefficient and ri{t) represents a thermal noise with zero mean 
and a correlator {ri(t)ri(t')) = 2'jkBT6{t — t'), T being the temperature of the medium. 
The force field F{x) = —d(j)/dx originates from the background potential 0(a;). We will 
consider the case when the potential 0(x) is a Brownian motion in x. This means the 
force F{x) is simply a Gaussian white noise in space with zero mean and a correlator 



The physics of this model is rather different in the two complimentary limits namely 
the overdamped and the underdamped case. In the overdamped limit where m = 0, 
this model reduces to the celebrated Sinai model|]l[. Since its introduction the Sinai 
model has been studied extensively over the past two decades[^. It is well known 
that the disorder drastically alters the diffusive behaviour of the system and the mean 
squared displacement at any nonzero temperature grows with time as (xl) ~ In'^(t) for 
large t, where the angle brackets here indicate the thermal average (over rf) and the 
overline indicates the disorder average (over F). 

In the underdamped limit (m > 0) this model represents, in the absence of thermal 
fluctuations (T = 0), the deterministic motion of a massive particle in a random 
medium. Recently there has been a lot of interest in understanding the dynamics of 
manifolds in a random mediumQ where one neglects the thermal fluctuations as a first 
approximation. The deterministic limit 17 = of Eq. (|l|) with m > then represents a 
toy model where one considers the motion of a zero dimensional manifold. For nonzero 
friction 7 > 0, the particle will eventually stop at some distance from the starting 
point since there is no thermal noise. Recently Jespersen and Fogedby studied the 
distribution of this stopping distance numerically and via simple scaling analysis Q. 
In the limit of vanishing friction 7 = 0, this model was earlier studied by Stepanow 
and Schultz0 who showed that the particle resumes normal diffusion in dimensions 
d > 1 provided it starts with a nonzero initial velocity vq. However, in one dimension 





{F{x)F{x')) =5{x-x'). 
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with 7 = 0, the particle does not diffuse but oscillates between two turning points as 
mentioned before. In this Letter we focus on this limiting case in one dimension and 
calculate the distribution of the oscillation period exactly. 

We set m = 1 for convenience and study simply the equation 

(2) 

where F{x) = —d(j)/dx is a white noise with zero mean and a correlator {F{x)F{x')) = 
6{x — x'). The potential (pix) = 00 + lo F{x')dx' is then a Brownian motion in x 
which starts at the value 0(0) = 0o at x = 0. Without loss of generality we set 
00 = 0. We assume that the particle starts at t = at x = with an initial velocity 
V > 0. The particle will move to the right until it reaches a turning point Xc where 
0(xc) = and then it will move back down to the point where its velocity will 
be —V (see Fig. 1). We define by T the time taken to arrive at the point Xc starting 
from X = 0. The particle will then move to the left till it similarly reaches the left 
turning point x'^, the time taken to do this is T'. The total period of oscillation is then 
given by Tosc = 2(T + T'). Since 0(x) is a Brownian motion in x, it follows from its 
Markov property that 0(a;) for a; > and x < are completely independent of each 
other. Thus T and T' are also independent and has the identical distribution which 
we compute below. The distribution of Tosc then follows simply as the convolution of 
the distribution of 2T with that of 2T'. 

It is easy to express T explicitly in terms of the potential 0(x), 

T= \ I (3) 
•^0 ^t;2 - 20(x) 

where we have used 0o = 0. Since 0(x) is random, the period T is also a random 
variable whose distribution we denote by P{T). Even though in the physical problem 
the turning time T is given by the formula in Eq. (|^) with 0o = 0, it is useful to consider 
the random variable T in Eq. (|^) as a function of an arbitrary 0o and compute the 
distribution P(T, 0o) that depends on 0o. Eventually the distribution of the physical 
time T is given by P{T) = P(T, 0). We define the Laplace transform 

Qifi, 0o) = / P(r, 0o) exp(-/ir) dT. (4) 

JO 

Substituting the explicit form of T from Eq. (^ in the Laplace transform one finds 



Q(/i, 0o) = E, 



<t>0 



exp (y—jj, J V{(l){x))dx^ 



(5) 



where E^^ denotes the expectation value over Brownian paths starting at 0o at time 
X = and 

V^(0(x)) = (6) 
^v^ - 20(x) 
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There is a nice probabilistic interpretation of (^(/i, 0o) in Eq- (H)- Consider a 
Brownian trajectory {(f){x) : < x < 00} starting at (pQ at x = 0. If the particle 
is killed with probability fiV{(t>{x))dx in the time interval [x,x + dx], then 0o) is 
simply the probability that the particle is not killed before reaching the level = 

In this interpretation Xc = where Ha is the first hitting time at level a, i.e., 

2 

Ha = infx>o{x : (f){x) = a}. 

To calculate Q{fi, 0o) we employ the standard backward Fokker Planck technique 
P, 1^ in the following way. One first evolves the Brownian motion over an infinitesimal 
time interval dx from the initial time x = 0. The Brownian path starts at (po at x = 0. 
In time dx it evolves to a new position 0o + (^(po- Thus the initial position of the 
trajectory that starts at time dx is 0o + d4>o- As a result of this infinitesimal change 
the function Q{fi, 0o) evolves as 



Q{fi, 4>o) = (1 - ^V{(f)o)dx) Ed4,a 00 + dcpo)] . 



(7) 



The Eq. (|^) follows from the fact that the particle survives in the interval [0, dx] with 
probability (1 — ^V{(j)Q)dx) and then restarts its trajectory from the point 0o + d'Po at 
time dx. Of course the variable dipQ is random and one needs to sum over all histories 
of evolution in the interval [0, dx] which is denoted by the expectation Ed^^ in Eq. (|^. 
Using the Ito prescription we average over dcpQ to first order in dx and then equating 
the terms of order dx yields the time independent Schrodinger equation 



1 rf2 



2d(i)l 



Q{fi, 0o) - fiV{(t>o)Q{fi, (po) 



0, 



The function Q satisfies the following boundary 

then 



where V{(p) is given by Eq. 

conditions: (i) Q{n,v'^/2) = 1. This follows from the fact that if (pQ = 

by definition the hitting time H^ = and subsequently Xc = in Eq. 

2 

Q{lJ., —00) = 0, as if 0(0) —00, the hitting time H^ — > 00 implying Xc ■ 

2 

Subsequently the integral in Eq. (^) diverges since V^(0(x)) > even when x 
and one obtains this boundary condition. 

I) we get 



00. 

■ 00 



Making the change of variables y 

dy'^ 



a/w^ — 200 in Eq. 
2fiyQ = 0. 



IdQ 
y dy 



(9) 



The general solution to this equation is 



Q{y) = AyKys 




Byh/z 




(10) 



where A and B are constants and Kv{z) and Iv{z) are modified Bessel functions of 
index v. As 0o —00, y — » 00. In this limit the boundary condition (ii) implies the 
constant 5 = since Iv{z) diverges for large argument. The boundary condition (i) 
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implies that as y ^ 0, Q{y) 1 which fixes the constant A = 2^/ V^^Vr (2/3)3^/3. 
Changing back to the 0o variable, we then have the exact expression for Q{fi, 0o) 



r,- - 

/3/i3 



0-^2/3 
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The distribution P{T, 0o) is then obtained by inverting the Laplace transform in Eq. 
(|Tl|) which fortunately can be done exactly 0. Setting 0o = we finally obtain 



P(T) 



exp 



2v^ 
'9T 



One can easily check that this distribution is normalized, 
able T = Ip- we find that the probability density p(r) - 
independent of v and is given by the simple expression 



p(r) 



exp — 



r 



(12) 

In terms of the scaled vari- 
= P{T)dT / dr of r becomes 



(13) 



which is valid for all t > 0. In particular for large r, p{t) ~ r~" decays as a power 
law with the exponent a = |. Thus none of the integer moments of the distribution 
exists and there are periods of all sizes extending upto infinity. The distribution P{T) 
has a unique maximum at T = = j^v^ indicating the existence of a most probable 
turning time Tm- 

The exponent a = 5/3 can also be derived by an intuitive scaling argument. /^From 
the energy conservation it follows that 0(xc) — ^(xo) = {dx^/ dtY /2 where Xc is the 
turning point, Xq is any other reference point and dx^/dt is the velocity at Xq. For 
large \xc — Xq\, clearly (p^Xc) — (^(xq) typically scales as ~ \l\xc — Xq\ since 0(x) is a 
Brownian motion. This indicates dxo/dt ~ \xc — xq\^^^ and hence \xc — xq] ~ f^^^ for 
large t. The distribution of Xc behaves as p{xc) ~ for large Xc since Xc is simply 

the first passage time to the level f^/2 of a Brownian walker starting at the origin. 
It then follows that the distribution of the turing time t ~ |xc — Xo|^^^ has a power 
law tail with a = 5/3. In fact this scaling argument can be generalized to an 
arbitrary random potential with a surface roughness exponent such that (pix) ~ x^ 
(x = 1/2 for the Brownian case). The distribution of the first passage time in general 
has a power law decay, p(xc) ~ x'^'^" where 9s is the spatial persistence exponent |0. 
Following the same argument as in the Brownian case, one finds that the distribution 
of turning time has a power law tail with an exponent a = 1 + 26^/(2 — x). For the 
Brownian case where x = 1/2 and Og = 1/2, one recovers a = 5/3. 

The distribution (|l^) has a rather remarkable link with the original (over damped) 
Sinai model in the interval [0, oo], in presence of an external force Fext = —fo in the 
interval [0, oo]. Here the partition function in the canonical ensemble is given by 



exp [—P{(j){x) + fox)] dx 



(14) 
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where (5 is the inverse temperature. For /o > the distribution p{Z) of Z exists, can 
be computed exactly [p!0[] and is given by 



Ml 

,3 



Hence, when fo = 13/3 the partition function Z (up to a constant multiphcative factor) 
for the statics of the Sinai model in the presence of a constant force — /o (towards the 
origin) has the same distribution as the turning time T in our model. 

We have also verified our exact results via the numerical integration of Eq. (^). We 
chose a discretization of the integral as follows 

T{^ + l) = T{i)+ (16) 

- 20(0 

0(i + l) = (p{i) + ai^Ki (17) 

Here the cxj are independent gaussian variables of mean zero and variance one. The 
Brownian limit is clearly approached when Ax 0. This system is integrated up 
till the point i where 0(i + 1) > ^ ^^^d the value of T{i) recorded. This is done for 
5 X 10^ realizations to construct the histogram of the T. The fact that T has a broad 
distribution implies that the number of integration steps required for a given sample 
can become very large. Therefore a temporal cutoff T* is placed on the values of T 
measured. This means we discard those samples where T exceeds a large fixed value 
T* . Thus the numerical histogram, once normalized, has the distribution 

^ P{T)e{T* - T) 
^ ' ^ P{S)dS ■ ^ ' 

In the limit when T* is large, P{T, T*) —>■ P{T). We first compute P(T, T*) analytically 



from Eq. (|l^) using the exact form of P{T) from Eq. (p!2|). We then compare this 
analytical result with the numerically obtained P{T,T*) in Fig. (0) using the scaled 
variable r = 9T/2v^ with the choice v = 1 and T* = 5000. The agreement is evidently 
excellent. 
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Figure 1: Particle in a Brownian potential with initial velocity v. The right/left 
turning points are shown as Xc and and x'^ respectively where the potential first becomes 
equal to 4> — 
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Figure 2: Histogram of r (circles) generated from 5 x 10^ samples, with initial velocity 
V = 1 and with temporal cutoff T* = 5000, the bin size is 0.1. Shown also is the 
analytical prediction. 
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